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Summary 

Most available knowledge on fungal arginine metabo- 
lism is derived from studies on Saccharomyces cer- 
evisiae, in wliich arginine catabolism is initiated by 
releasing urea via the arginase reaction. Orthologues 
of the S. cerevisiae genes encoding the first three 
enzymes in the arginase pathway were cloned from 
Kluyveromyces lactis and shown to functionally com- 
plement the corresponding deletion in S. cerevisiae. 
Surprisingly, deletion of the single K. lactis arginase 
gene KICAR1 did not completely abolish growth on 
arginine as nitrogen source. Growth rate of the dele- 
tion mutant strongly increased during serial transfer 
in shake-flask cultures. A combination of RNAseq- 
based transcriptome analysis and ^^C-^^N-based flux 
analysis was used to elucidate the arginase- 
independent pathway. Isotopic ^^C^^N-enrichment in 
y^aminobutyrate revealed succinate as the entry point 
in the TCA cycle of the alternative pathway. Transcript 
analysis combined with enzyme activity measure- 
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ments indicated increased expression in the KIcarIA 
mutant of a guanidinobutyrase (EC. 3.5.3.7), a key 
enzyme in a new pathway for arginine degradation. 
Expression of the K. lactis KLLA0F27995g (renamed 
KIGBU1) encoding guanidinobutyrase enabled S. cer- 
evisiae to use guanidinobutyrate as sole nitrogen 
source and its deletion in K. lactis almost completely 
abolish growth on this nitrogen source. Phylogenetic 
analysis suggests that this enzyme activity is wide- 
spread in fungi. 

Introduction 

Arginine metabolism lias been subject of intensive bio- 
cliemical studies. After discovery of the urea cycle for use 
of arginine as a nitrogen source (Krebs and Henseleit, 
1932; Krebs, 1973), attention focused on its role as a 
precursor for the synthesis of polyamine and the signal- 
ling compounds y-aminobutyrate (GABA) and nitric oxide 
(Knowles and IVIoncada, 1994; Pitkanen etal., 2001). The 
most widely distributed pathway for arginine degradation 
that occurs across all three kingdoms (Abdelal, 1979) is 
initiated by arginase (EC 3.5.3.1), a ureohydrolase that 
converts arginine to ornithine and urea. Its active site, 
which contains several Mn^+ binding sites, is also con- 
served in other ureohydrolases such as agmatinase (EC 
3.5.3.11), formiminoglutamase (EC 3.5.3.8), proclavami- 
nate amidinohydrolase (EC 3.5.3.22), guanidinobutyrase 
(EC 3.5.3.7) and guanidinopropionase (EC 3.5.3.17) 
(Ouzounis and Kyrpides, 1994). Genes encoding these 
enzymes are assumed to have emerged early in evolution 
(Hartman, 1975) and have been used as markers in phy- 
logenetic studies (Ouzounis and Kyrpides, 1994; 
Sekowska etal., 2000). 

In eukaryotes, only two types of ureohydrolase have 
hitherto been described. In addition to arginase, higher 
eukaryotes express agmatinase (Coleman etal., 2004), 
which participates in an alternative pathway for arginine 
catabolism (Fig. 1). In this pathway, arginine is first decar- 
boxylated to agmatine, which is converted to putrescine 
and urea by agmatinase. Putrescine can then be con- 
verted either to GABA or to the polyamines spermine and 
spermidine (Pegg, 2009). The rapidly Increasing number 



© 2014 The Authors. Molecular Microbiology published by John Wiley & Sons Ltd. 

This is an open access article under the terms of the Creative Commons Attribution-NonCommercial-NoDerivs License, which permits use 
and distribution in any medium, provided the original work is properly cited, the use is non-commercial and no modifications or adaptations 
are made. 



370 G. Romagnoli et al. 



L-argino succinate 

H 



EC.6.3.4.5 



OH 

y-aminobutyrate 
I EC.1.2.1.3 

y-aminobutanal 

EC.1.4.3.10 




EC4.1.1.17 



EC.2.6.1.13 

o 



HO 

L-glutamate-Y-semialdehyde 



spontaneous 

„ ? EC.1.5.1.2 H o 

^ {^-''^- 

L-A^-pyrroline-S-carboxylate Proline 
^ EC.1. 5.99.8 I 



OH 

NHj 

Glutamate 



EC.1.5.1.12 



O 

Q. 



3 

CD 

o 
3 



Fig. 1. Overview of the key reactions in eul<aryotic arginine metabolism. Thicl< lines indicate ureohydrolase reactions. EC 1.2.1.3: aldehyde 
dehydrogenase:, EC 3.5.3.11: agmatinase, EC 4.1.1.19: arginine decarboxylase, EC 4.3.2.1: argininosuccinate lyase, EC 6.3.4.5: 
argininosuccinate synthase, EC .1.14.13.39: nitric oxide synthase, EC 1.4.3.10: putrescine oxidase, EC 2.1.3.3: ornithine 
carbamoyltransferase, EC 3.5.3.1: arginase, EC 4.1.1.17: ornithine decarboxylase, EC 2.6.1.13: ornithine aminotransferase, EC 1.5.1.2: 
pyrroline-5-carboxylate reductase, EC 1.5.99.8: proline dehydrogenase, EC 1.5.1.12: L-A'-pyrrollne-5-carboxylate dehydrogenase. 



of wliole-genome sequences has enabled the putative 
identification of arginase and agmatinase genes in many 
eukaryotes. However, since such annotation is based on 
sequence homology only, it does not enable definitive 
conclusions on the catalytic function of the encoded 
proteins. 

Much of the knowledge on fungal arginine metabolism 
is based on studies with the model organism Saccharo- 
myces cerevisiae. In S. cerevisiae, arginine is transported 
into the cell and subsequently hydrolysed by arginase 
(Carl) to yield ornithine and urea (Cooper etal., 1992; 
Sumrada and Cooper, 1 992; Shima et al., 2003). An ATP- 
dependent amidolyase (Dur1,2) then converts urea into 
ammonia and carbon dioxide. Ornithine is further con- 
verted by an ornithine-specific transaminase (Car2) into 
glutamate-y-semialdehyde (GSA), which spontaneously 
forms 1-pyrroline-5-carboxylate (P5C) (Martin etal., 



2003). Due to subcellular compartmentation, S. cerevi- 
siae is unable to convert cytosolic P5C directly to gluta- 
mate (Davis, 1986). Instead, P5C is reduced to proline 
using pyroline-5-carboxylate reductase (Pro3). Proline is 
then transported into the mitochondria (Brandriss and 
Falvey, 1992), converted back to P5C by an oxidase 
(Puti) and, finally converted to glutamate by mitochon- 
drial P5C dehydrogenase (Put2) (Davis, 1986). Since only 
very few physiological studies have been conducted on 
arginine metabolism in non- Sacchammyces yeasts, it is 
unknown whether the arginase pathway, which is essen- 
tial for growth of S. cerevisiae on arginine as sole nitrogen 
source (Bossinger and Cooper, 1977), is the only fungal 
pathway for arginine catabolism. 

Sacciiaromyces cerevisiae and Kluyveromyces lactis 
both belong to the Saccharomycetaceae family. These 
two related yeasts are considered to have genetically 
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separated before the whole-genome duplication (WGD) 
event that reshaped the genome of S. cerevisiae, further- 
more K. lactis is regarded as resembling a pre-WGD 
ancestor of S. cerevisiae (Dujon, 2010). While many 
studies have been conducted on the differences in sugar 
metabolism between these two species, the differences in 
amino acid metabolism have not been studied in detail. 
Nonetheless, the complete genome sequence of K. lactis 
revealed many putative orthologues of S. cerevisiae 
genes involved in arginine metabolism (Dujon etal., 2004; 
Souciet etal., 2009; Dias etal., 2012). 

The goal of the present study was to characterize argi- 
nine catabolism in K. lactis. To this end, CAR1, CAR2 and 
PR03 orthologues in K. lactis were identified and func- 
tionally analysed by deletion, expression in S. cerevisiae 
and enzyme activity assays. Since deletion of the argin- 
ase gene in K. lactis was found not to completely abolish 
growth on arginine as a sole nitrogen source, the alterna- 
tive pathway for arginine catabolism operating in this 
yeast was studied by a combination of transcriptome 
analysis, ^^C and ^^N isotope-based flux analysis and 
enzyme activity assays in cell extracts. 

Results 

Arginine catabolism in K. lactis occurs primarily via the 
arginase pathway 

To investigate arginine metabolism in K. lactis, Blastp 
searches were performed to identify orthologues of the S. 
cerevisiae L-arginase (Carl), (Brandriss and l\/lagasanil<. 



1980; Sumrada and Cooper, 1992), ornithine transami- 
nase (Car2) (Dubois etal., 1978; Parl< etal., 1999), 
and pyrroline-5-carboxylate reductase (Pro3) (Brandriss 
and Falvey, 1992) in the K. lactis genome. The pre- 
dicted proteins KLLA0F00704p, KLLA0B05247p and 
KLLA0D10736p showed 65, 71 and 71% identity at the 
amino acid level to the S. cerevisiae proteins and were 
therefore tentatively named KICAR1, KICAR2 and 
KIPR03 respectively. Single-deletion strains for each of 
the three S. cerevisiae genes and for their putative K. 
lactis orthologues were constructed. Moreover, KICAR1, 
KICAR2 and KIPR03 were cloned in a S. cerevisiae 
expression vector and transformed to the corresponding 
S. cerevisiae deletion strains to test heterologous 
complementation. 

To verify whether deletion resulted in elimination of the 
corresponding enzyme activity, enzyme activity assays 
were performed on the mutant strains. After an initial 
growth phase on ammonium as sole nitrogen source [with 
the exception of the cells harvested for measurement of 
the P5C reductase activity which required addition of 
proline due to their proline auxotrophy (Tomenchol< and 
Brandriss, 1987)], cells were washed and incubated for 
2 h on arginine or ornithine as sole nitrogen source to 
induce arginine-related enzymes. Deletion of CAR1, 
CAR2 and PROS orthologues led to complete elimination 
of the corresponding enzyme activity in both K. lactis and 
S. cerevisiae (Table 1). Moreover, S. cerevisiae deletion 
strains expressing KICAR1, KICAR2 and KIPR03 ortho- 
logues showed a high activity for each of the correspond- 
ing enzymes (Table 1). 



Table 1 . Arginase, ornithine transaminase (OTA) and P5C-reductase enzyme activities measured in cell extracts of the yeast strains K. lactis 
GG1632 {KICAR1, KICAR2, KIPR03}, S. cerevisiae IMZ312 {CAR1, CAR2, PR03), K. lactis IMK432 {KIcarIA), S. cerevisiae IMK435 {car1A), 
S. cerevisiae IMZ310 {car1A KICARiT), K. lactis IMK433 {Klcar2A), S. cerevisiae IMK436 {car2A), S. cerevisiae \MZ3^^ (car2A KICAR2T), K. lactis 
IIV1K434 {Klpro3A), S. cerevisiae IMK445 (pro3A) and S. cerevisiae IMK342 (proSA KIPROST). 



Strain Activity (jimol min ' mg protein^') 



Arginase 


GG1632 {KICAR1) 


K. lactis 


1.45 ±0.3 




IMK432 (KlcarU) 


K. lactis 


B.D. 




IMZ312 (CAR1) 


S. cerevisiae 


3.7310.8 




IMK435 \car1A) 


S. cerevisiae 


B.D. 




IMZ310 (carl A KICAR1'[) 


S. cerevisiae 


12.59 ±3.0 


Ornithine transaminase 


GG1632 (KICAR2) 


K. lactis 


30.18 ±2.1 




IMK433 (Klcar2A) 


K. lactis 


B.D. 




IMZ312 (CAR2) 


S. cerevisiae 


421 .93 ±3.4 




IMK436 (car2A) 


S. cerevisiae 


B.D. 




IMZ311 (car2A KICAR27) 


S. cerevisiae 


606.28 ± 37.8 


P5C reductase 


GG1632 (KIPR03) 


K. lactis 


0.58 ±0.1 




IMK434 (Klpro3A) 


K. lactis 


B.D. 




IMZ312 (PROS) 


S. cerevisiae 


0.096 ± 0.0 




IMK445 (pro3A) 


S. cerevisiae 


B.D. 




IMZ342 (pro3A KIPR03T) 


S. cerevisiae 


25.71 ±2.0 



Cells were harvested after incubation with arginine, ornithine or proline. Data represent the average ± mean deviation of independent biological 
duplicates. B.D.: below detection limit of the enzyme assays, estimated at 0.07, 0.03 and 0.02 i^mol min ' mg protein ' for the arginase, OTA or 
and P5G-reductase measurement respectively. 
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Fig. 2. Phenotypic analysis of K. lactis GG1632 {KICAR1, KICAR2, KIPR03), IMK432 {KIcarIA), IMK433 {Klcar2A), IIVIK434 {Klpro3A) and S. 
cerevisiae IMZ312 (CAR1, CAR2, PROS), IMZ310 {carl A KICARit), IMZ311 {car2A KICAR21), IMK342 (proSA KIPR03\), IMK435 {carl A), 
IIV1K436 {car2A), IMK445 {pro3A) for growth on various nitrogen sources. The strains were plated on synthetic medium with glucose and either 
ammonium, arginine, ornithine or proline. The plates were incubated at 30°C and scored after 72 h. 



Phenotypes of the deletion strains and of the heterolo- 
gously complemented S. cerevisiae strains were first 
analysed on spot plates containing arginine, ornithine, 
proline or ammonia as a sole nitrogen source (Fig. 2). 
Consistent with earlier studies (Brandriss and Magasanik, 
1980; Brandriss and Falvey, 1992), the S. cerevisiae 
strains IMK435 {car1A), IMK436 (car2A) and IMK445 
{pro3A) were unable to grow on arginine as a sole nitro- 
gen source. Furthermore, growth on ornithine was abol- 
ished for IMK436 (car2A) and also for IMK445 {pro3A), 
consistent with the proline auxotrophy of pro3 mutants 
(Tomenchok and Brandriss, 1987). On spot plates, the 
same pattern was observed for K. lactis strains IMK432 
(KIcarIA) and IMK434 (Klpro3A) (Fig. 2). In contrast, the 
strain IMK433 (KIcarZA) exhibited a very slow but notice- 
able growth on ornithine plates. While unexpected this 
results might reflect a difference of specificity and higher 
substrate promiscuity of the K. lactis transaminases. Fur- 
thermore KICAR1, KICAR2 or KIPR03 \n the S. cerevisiae 
deletion mutants showed similar growth characteristics to 
the prototrophic empty-vector reference strain S. cerevi- 
siae ll\/IZ312, indicating that the K. lactis orthologues 
fully complemented the corresponding null mutations in 
S. cerevisiae. 

Growth of arginase-negative K. lactis on arginine as 
sole nitrogen source 

In contrast to the apparent absence of arginase activity 
and growth on solid medium, the K. lactis strain IMK432 



{KIcarIA) was still able to grow in liquid synthetic medium 
containing arginine as sole nitrogen. Growth in these cul- 
tures was characterized by an extended lag phase (esti- 
mated specific growth rate of 0.03 h"^), and it took about 6 
days to reach stationary phase (Fig. 3). During four con- 
secutive serial transfers the specific growth rate gradually 
increased up to 0.1 h"^ and the fully grown culture 
reached a final ODeoo of 30 in stationary phase. No growth 
was observed in serial transfer experiments with S. cer- 
evisiae IMK435 {car1A). To investigate whether the 
observed increase in growth rate was due to adaptation of 
the culture or was genetically stable, cells were trans- 
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Fig. 3. Adaptive growth profile of the strain IMK432 {KIcarIA Arg+) 
in sequential shake flask cultivations. The strain was grown in four 
successive shake flaks with chemically defined medium with 
glucose and arginine as nitrogen source. A single colony isolate 
from the last shake flask was selected, and named IMS0367. 
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Fig. 4. Growth kinetics of K. lactis GG1632 (left) and IIV1S0367 KIcarIA (right). 

A. Optical density at 660 nm (blue) and extracellular concentration of arginine (green), glucose (red), urea (purple) and ammonia (orange) 
were followed in aerobic batch cultures on glucose synthetic medium with arginine as nitrogen source. 

B and C. The total Carbon (B) and Nitrogen (C) pools were determined by HPLC and off-gas measurements taken during the exponential 
growth phase (15 h for GG1632 and 52 h for II\/1S0367) based on standard biomass composition estimation {Supplemental information) and 
are shown as a percentage of the initial pool. 



ferred to a shake flask containing glucose synthetic 
medium with ammonium sulphate as the nitrogen source 
and grown for approximately 20 generations. After rein- 
troduction into arginine medium, this culture did not show 
any delay or decrease in growth rate compared to the 
fourth transfer on arginine, suggesting that a mutation 
contributed to the fast growth of the arginine-adapted 
cultures. For further characterization, a single-colony 
isolate was obtained from the final transfer on arginine. 
The resulting K. lactis strain was named IMS0367 
(KIcarIA Arg+) and exhibited a specific growth rate of 
0.09 h~\ Deletion of the KICAR1 gene in this strain was 
confirmed by diagnostic PGR. 

Physiological and transcriptome analysis of an 
arginase-negative K. lactis strain 

To further investigate arginine metabolism in the evolved 
arginase-negative K. lactis strain, strains GG1632 



(KlkuSOA KICAR1 reference strain) and IMS0367 
{KIcarIA Arg+) were grown in aerobic bioreactor batch 
cultures on chemically defined medium with glucose and 
arginine as sole nitrogen source (Fig. 4). The specific 
growth rate of the arginase-negative strain was around 
fourfold lower than that of the reference strain (Table 2). 
The lower apparent biomass yield on glucose of the 
mutant strain (0.41 g-g"^ compared to 0.60 g-g"'' for the 
reference strain) was at least partially due to the fact 
that, in contrast to the reference strain; it did not com- 
pletely consume the available arginine (Table 2). Moreo- 
ver, no extracellular urea and/or ammonia were found 
during growth of IMS0367. For both strains, approxi- 
mately 98% of the carbon initially present in the form of 
arginine and glucose could be recovered at the end of 
cultivation as CO2, and biomass for GG1632 and as 
CO2, biomass and residual arginine, in case of IMS0367 
(Fig. 4). Similar observations could be drawn from the 
nitrogen balance indicating that, also in the evolved 



© 2014 The Authors. Molecular Microbiology published by John Wiley & Sons Ltd., Molecular Microbiology, 93, 369-389 



374 G. Romagnoli et al. ■ 

Table 2. Physiological data of K. lactis strains GG1632 and IMS0367 {KIcarIA Arg*) grown in batch cultures on synthetic glucose medium with 
arginine as sole nitrogen source. 



GG1632 IMS0367 



Growth rate (h-') 


0.38 ±0.01 


0.087 ±0.001 


Ysx.giuc (gdw-g ') 


0.60 ± 0.04 


0.41 ±0.01 


Arginine consumed (mM): 


19.8 ±0.0 


8.2 ±0.1 


Dry Weight (g M): 


10.8 ±0.5 


7.2 ±0.36 


Approx. qs.giuc (g • h '-gdw"') 


0.57 ± 0.07 


0.22 ± 0.07 


Approx qs.ara (g-h '-gdw ') 


0.19 ±0.02 


0.012 ±0.003 


Arginase activity ( |imol min^' mg protein^') 


0.516 ±0.15 


B.D. 



The batch cultures were sampled in mid-exponential phase, cell extracts were prepared and arginase activity was measured. Data represent the 
average ± mean deviation of independent biological duplicate cultures and enzymatic measurements. B.D. denotes below the detection limit, which 
was estimated at 0.07 )j.mol min^' mg protein 



mutant strain, all four nitrogen atoms in arginine were 
utilized (Fig. 4). 

RNA sequencing of samples taken during the exponen- 
tial phase of growth on glucose-arginine media of the 
reference strain and the evolved arginase-negative strain 
revealed 430 genes whose transcripts levels differed by 
more than twofold between the two strains. Functional 
category enrichment analysis (Kanehisa etal., 2008; 



Knijnenburg etal., 2009; Hunter etal., 2012) using KEGG 
and InterPro database annotations identified significant 
over-representation of genes involved in arginine and 
proline metabolism, amino acid transporters and alanine, 
aspartate and glutamate metabolism (Table 3). The tran- 
scriptional data mimicked the cellular response to nitro- 
gen limitation and/or growth on a non-preferred nitrogen 
source. This is illustrated by the upregulation of the 



Table 3. KEGG and INTERPRO categories overrepresented in the set of genes that showed significantly different (P-value < 0.005) transcript 
levels in K. lactis IMS0367 {KIcarIA Arg*) relative to the reference strain K. lactis GG1632 1(FC > 121, FDR 1%). 



Arginine and proline metabolism (KEGG, P-value: 5E-7, 15 of 65 genes) 



FC 



KLLA0F27995g 


Conserved hypothetical protein/putative arginase 


43.4 


KLLA0B08800g 


Similar to S. cerevisiae YDR242W AMD2 putative amidase 


14.7 


KLLA0F04235g 


Conserved hypothetical protein/putative arginase 


13.5 


KLLA0E08119g 


Highly similar to S. cerevisiae YBR208C DUR1 2 urea amidolyase 


6.1 


KLLA0F14366g 


Similar to S. cerevisiae ARG3 ornithine carbamoyltransferase 


4.1 


KLLA0F00594g 


Highly similar to S. cerevisiae GDH3 NADP(+)-dependent glutamate dehydrogenase 


3.9 


KLLA0B05247g 


Highly similar to S. cerevisiae CAR2 L-ornithine transaminase (OTAse) 


2.9 


KLLA0E21077g 


Highly similar to S. cerevisiae ARG8 acetylornithine aminotransferase 


2.8 


KLLA0F26268g 


Similar to S. cerevisiae mitochondrial ornithine acetyltransferase catalyses 


2.2 


KLLA0C07997g 


Highly similar to S. cerevisiae ARG5 enzyme 


2.2 


KLLA0A08492g 


Conserved hypothetical protein 


2.0 


KLLA0D16962g 


Similar to S. cerevisiae PUT1 proline oxidase 


-2.5 


KLLA0C1 0549g 


Highly similar to S. cerevisiae PUT2 delta-1 -pyrroline-5-carboxylate dehydrogenase 


-7.3 


Amino acid transporters (INTERPRO, P-value: 0.003, 7 of 323 genes) 


KLLA0B14685g 


Weakly similar to S. cerevisiae PUT4 proline-specific permease 


6.4 


KLLA0C02343g 


Similar to similar to S. cerevisiae CAN1 plasma membrane arginine permease 


6.3 


KLLA0A06930g 


Similar to S. cerevisiae GAP1 general amino acid permease 


4.0 


KLLA0F01012g 


Similar to S. cerevisiae GNP1 high-affinity glutamine permease 


2.8 


KLLA0A10813g 


Similar to S. cerevisiae TAT2 tryptophan permease high affinity 


2.2 


KLLA0C01606g 


Similar to S. cerevisiae GNP1 high-affinity glutamine permease 


-6.9 


KLLA0A06886g 


Similar to S. cerevisiae GAP1 general amino acid permease 


-15.7 



Alanine, aspartate and glutamate metabolism (KEGG, P-value: 0.005, 6 of 65 genes) 



KLLA0F20548g 


Similar to S. cerevisiae UGA1 gamma-aminobutyrate (GABA) transaminase 


4.7 


KLLA0F24882g 


Similar to S. cerevisiae CPA1 small subunit of carbamoyl phosphate synthetase 


4.1 


KLLA0F00594g 


Similar to S. cerevisiae GDH3 NADP(+)-dependent glutamate dehydrogenase 


3.9 


KLLA0F13640g 


Similar to S. cerevisiae mitochondrial aspartate aminotransferase 


3.4 


KLLA0D12980g 


Highly similar to S. cerevisiae ASP1 cytosolic L-asparaginase 


-2.2 


KLLA0C1 0549g 


Highly similar to S. cerevisiae PUT2 delta-1 -pyrroline-5-carboxylate dehydrogenase 


-7.3 



Both strains were grown on batch cultures with glucose and arginine as sole N-source. P-values were calculated using Fischer exact statistics. 
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Table 4. Intracellular concentrations, and corresponding final 


"C'=N enrichments (%) 


of arginine catabolism intermediates measured in K. lactis 


IMS0367 (KIcarIA Arg*) and K. lactis GG1632 (reference strain, KLCAR1 Arg+). 












Labelled isotopologue species 






Intracellular cone. (|j,mol gdw 




(final overall "C"N enrichment) 


Metabolite 


GG1632 


IMS0367 


GG1632 


IMS0367 


Arginine 


370 ± 80 


660 ± 72 


38% 


26% 


Proline 


7.24 ± 1.92 


0.06 ± 0.02 


83% 


9% 


Glutamate 


158 ±47 


48± 17 


86% 


44% 


a-Ketoglutarate 


2.9710.81 


2.28 ± 0.62 


55% 


3% 


Succinate 


3.01 ±1.2 


6.61 ±0.93 


49% 


5% 


GABA 


0.01 ±0.01 


0.41 ± 0.33 


3% 


42% 


Putrescine 


2.56 ± 0.89 


0.03 ±0.02 


27% 


B.D. 



Cellular extracts were prepared from batch cultures grown on glucose synthetic medium with fully ("C'^N) labelled arginine as sole nitrogen source. 
Data represent the average ± mean deviation of independent biological duplicate cultures. B.D. denotes below the detection limit, N.D. denotes 
not determined. 



K. lactis orthologues of several genes that, in S. cerevi- 
s/ae are under nitrogen catabolite repression, amino acid 
transporters, and amino acid catabolism-related genes 
(Boer etal., 2005; 2007). Interestingly, tlie K. lactis 
orthologues of S. cerevisiae PUT1 and PUT2 genes 
involved in proline degradation (KLLA0D1692g and 
KLLA0C10549g), were significantly downregulated in 
IMS0367. 

Among the genes that were transcriptionally upregu- 
lated in K. lactis IMS0367, two open reading frames 
(KLLA0F27995g and KLLA0F04235g) shared sequence 
similarity with ureohydrolase encoding genes. As no argi- 
nase activity was detected in K. lactis IMS0367 (Table 2), 
it was hypothesized that the proteins encoded by these 
two genes might catalyse a similar reaction but with a 
different substrate. Several genes related to GABA 
metabolism were also upregulated in K. lactis IMS0367. In 
addition to the upregulation of KLLA0E14675g, which is 
similar to the S. cerevisiae PUT4 gene encoding the 
GABA-accepting permease, the transcription of a putative 
y-aminobutyrate transaminase, which in S. cerevisiae is 
responsible for the initial step in GABA catabolism 
(Vissers etai., 1989; Bach etal, 2009), was 4.7-fold 
upregulated in K. lactis IMS0367 (Table 3). 

Carbon and nitrogen isotopologue labelling using 
^^C^^N L-arglnlne 

To investigate whether the alternative route of arginine 
degradation was active in K. lactis IMS0367 {KIcarIA 
Arg+), a ^^C^^N labelling experiment was performed in 
shake flask cultures with fully labelled '^C6^^N4 arginine 
(M+'°) as a sole nitrogen source. Intracellular metabolite 
concentrations and ^^G^^N enrichment levels for each 
metabolite were determined over the course of 64 min of 
cultivation for strain GG1632 cultivation and 300 min for 
strain IMS0367 (Figs S1 and S2). The already high intra- 



cellular arginine concentration in strain GG1632 [370 
l^mol g"^ (biomass dry weight), corresponding to an 
intracellular concentration of 185 mlVI assuming an intra- 
cellular volume of c. 2 ml g"^] was found to be even 
higher in IMS0367 at 660 |xmol g"^ (biomass dry weight). 
Combined with the effect of vacuolar storage (Davis, 
1986) this resulted in an overall arginine enrichment that 
only reached 38% and 26% for strains GG1632 and 
IMS0367 after 64 and 300 min respectively (Table 4 and 
Table SI). 

In the reference strain K. lactis GG1632 (KiCARI), the 
enrichment found for fully labelled ^^C'^N isotopologue 
species of '^Cb^^N proline, ^^Cb'^N glutamate and ^^C5^^N2 
glutamine were in agreement with the presence of the 
arginase pathway described for S. cerevisiae as they 
could only derive directly from the ^^C6^^N4 labelled argi- 
nine (Fig. 5). In contrast, the enrichment of the carbon 
labelled isotopologue species in strain IMS0367 (KIcarIA 
Arg+) were found to reach a maximum level of maximally 
4%, thereby confirming that the arginase pathway was 
inactive. A substantial enrichment of only the ^^N gluta- 
mate isotopologue species confirmed that the cells were 
able to mobilize the nitrogen contained in ^^C^^N arginine. 
While for K. lactis GG1632 a distinct enrichment in fully 
labelled ^^Cs a-ketoglutarate was observed, the absence 
of such an isotopologue species, combined with the iden- 
tification of fully labelled succinate, fumarate and malate, 
indicated that, in K. lactis IMS0367, carbon originating 
from arginine entered the TCA cycle at a different point. 
Isotopologue measurements of GABA subsequently 
revealed that, unlike the situation in K. lactis GG1632, a 
substantial increase in fully labelled '^Cs^^N GABA was 
found in strain IMS0367. In eukaryotes, GABA is an inter- 
mediate of an arginine degradation pathway that also 
involves the polyamine putrescine. In K. lactis GG1632, 
enrichment of fully ^^C4'^N2 labelled putrescine was 
observed, most likely originating from decarboxylation of 
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Fig. 5. Enrichment of intracellular metabolites and flux distribution of K. lactis GG1632 {KICAR1) and K. lactis IMS0367 {KIcarIA Arg*) grown 
on "C6'^N4 arginine. Enrichment profiles are shown next to the corresponding compounds in plots with light blue background for strain 
GG1632 and green for the IMS0367 strains. Red lines refer to the species only. The fluxes are indicated in boxes (expressed in |imol 
gdw"' min fluxes indicated in blue correspond to strain GG1632 (KICAR1 Arg*) and those indicated in green to IIVIS0367 {KIcarIA Arg+). 
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Table 5. Guanidinobutyrase activities measured in cell extracts of K. lactis and S. cerevisiae strains GG1632, IMS0367, IME215, IME216 and 
CEN.PK113-7 D grown in batch cultures with arginine as sole nitrogen source. 

Activity 

Strain Description (limol min ' mg protein ') 



K. lactis GG1632 Prototrophic reference KICAR1 KIGBU1 B.D. 

K. lactis IMS0367 KIcarIA Arg* 0.14 ± 0.004 

S. cerevisiae IME215 MATa ura3-52 pUDE264 {TDH3p,-KIGBU1-CYC1,e, URA3) 0.17 ± 0.006 

S. cerevisiae IME216 MATa ura3-52 pUDE265 {TDH3p,- KLLA0F04235g-CyCyter URA3) B.D. 

S. cereWs/ae CEN.PK113-7D Prototrophic reference B.D. 



The K. lactis strains GG1632 and IMS0367, as well as the S. cerevisiae strains IME216 and IME215 were pre-grown in synthetic medium with 
glucose and ammonium as sole nitrogen source. B.D. denotes below detection limit, which was estimated at 0.005 |jmol min ' mg protein Data 
represent the average ± mean deviation of independent biological duplicate cultures. 



^^C5'^N2 ornithine. Intracellular putrescine levels in K. 
lactis IMS0367 were 10-fold lower than in K. lactis 
GG1632. Moreover, we were unable to quantify any of the 
isotopologue species of putrescine, as the signal was 
below the detection limit. For arginine, a slight increase in 
the ^^Ce^^Ns (IVI*^) isotopologue species was observed in 
both GG1632 (3%) and IIVIS0367 (6%), probably caused 
by reversible arginine transamination. 

To determine whether arginine catabolism in K. lactis 
IMS0367 indeed involved GABA as an intermediate, a 
metabolic flux model of the hypothetical pathway was 
constructed (Fig. 5). In this model, an alternative catabolic 
route was evaluated, in which arginine is first transami- 
nated, yielding ketoarginine and glutamate. Ketoarginine 
is then transformed into GABA via a decarboxylation and 
subsequent hydrolysis, releasing CO2 and urea respec- 
tively. The model was not solely based on a ^^C-calculated 
flux analysis, because the co-consumption of unlabelled 
glucose prevented accurate and unequivocal quantifica- 
tion of all intracellular metabolic fluxes. Consequently, the 
model consisted of two separate parts, in which metabolic 
fluxes were quantified based on data obtained from nitro- 
gen (Fig. S3 and Table S2) and carbon labelling (Fig. S4 
and Table S3) of arginine respectively. To account for the 
impact of vacuolar storage of arginine on label enrichment 
patterns, separate cytosolic and vacuolar pools were 
introduced in the model. The observed increases in flux 
through TCA cycle reactions were in agreement to the 
previously observed arginine and glucose uptake rates 
during the batch cultivations. For K. lactis GG1632, the 
metabolic flux analysis demonstrated that nitrogen origi- 
nating from arginine was not solely incorporated into 
biomass. Instead, 'excess' nitrogen was excreted in the 
form of ammonium, as observed experimentally during 
batch cultivation (Fig. 4). The flux analysis of K. lactis 
IMS0367 cultivation showed a different pattern. The 
model was only able to accurately describe the observed 
enrichment profiles when it was assumed that all nitrogen 
contained in the arginine molecule was mobilized in form 
of urea or glutamate. Flux calculations focusing on the 



fate of the carbon skeleton of arginine revealed that the 
rate of enrichment of the TCA cycle intermediates of the 
arginase-negative strain K. lactis IMS0367 strain was fully 
compatible with the degradation of GABA into succinate 
(Fig. 5, Figs S2 and S4). 

Identification of a functional guanidinobutyrase in 
K. lactis 

Transcriptome data and '^C^^N flux analysis identified 
GABA as a key intermediate in arginine catabolism in the 
arginase-negative strain K. lactis II\/1S0367 (KIcarIA Arg*). 
Combined with increased transcript levels of two 
putative ureohydrolase genes (KLLA0F27995g and 
KLLA0F04235g) in this strain, this strongly suggested 
involvement of a ureohydrolase in the alternative arginine 
catabolic pathway. To test this hypothesis, we measured 
the ureohydrolase activity with both agmatine and guani- 
dinobutyrate as substrate in cell extracts of the K. lactis 
strains IMS0367 and GG1632, grown on glucose with 
arginine as sole nitrogen source (Table 5). No activity with 
either substrate was measured in the reference strain 
GG1632, while the arginase-negative strain IMS0367 
exhibited a high activity exclusively with guanidinobu- 
tyrate (Table 5). Additionally, in line with these enzyme 
activity measurements, shake flask experiments revealed 
that IMS0367 grew on media containing guanidinobu- 
tyrate as sole nitrogen source. The K. lactis reference 
strain GG1632 was also able to use guanidinobutyrate as 
nitrogen source, indicating that the expression of the 
guanidinobutyrase activity was tightly controlled and at 
least induced in presence of guanidinobutyrate in KICAR1 
strains. In contrast, the prototrophic S. cerevisiae strain 
CEN.PK113-7 D was unable to use guanidinobutyrate as 
nitrogen source. 

To identify which of the two putative ureohydrolase 
genes in K. lactis contributed to the measured guanidi- 
nobutyrase activity, KLLA0F27995g and KLLA0F04235g 
were cloned under the control of the strong constitutive 
TDH3 promoter {TDH3pr) in an expression vector and 
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Fig. 6. Typical growth profile of K. lactis and S. cerevisiae strains 
on guanidinobutyrate: The K. lactis strains GG1632 (KICAR1) (•), 
II\/1S0367 {KIcarIA) (T), IMK560 {Kicar1AKigbu1A) (■) and the S. 
cerevisiae strain IME215 {TDH3pr::KiGBU1::CYC1,er) (♦) were 
grown in shake flasks on glucose synthetic medium containing 
guanidinobutyrate as sole nitrogen source. 



transformed to S. cerevisiae. Cell extracts of the result- 
ing strains did not exhibit any agmatinase activity 
but only cell extracts of S. cerevisiae IME215 
(rD/-/3pr::KLLA0F27995g) showed a substantial activity 
with guanidinobutyrate, cell extracts of S. cerevisiae 
IME216 {TDH3prKLLA0F04235g-CYC1,er) did not exhibit 
any measurable activity (Table 5). In agreement with these 
results, S. cerevisiae IME215 (7'DH3pr::KLLA0F27995g) 
was able to grow on guanidinobutyrate as sole nitrogen 
source (Fig. 6), resulting in a specific growth rate of 
approximately 0.1 h"^ in chemically defined glucose 
medium. Therefore, the gene KLLA0F27995g was 
renamed KiGBUL 

To further confirm the function of KLLA0F27995g/ 
KiGBUI and its role in an arginase-independent arginine 
catabolism pathway, KLLA0F27995g was deleted in 
IMS0367. While the K. lactis reference strain GG1 632 and 
the KicarIA Arg+ strain IMS0367 both grew at a specific 
growth rate of 0.45 h"^ on medium containing guanidi- 
nobutyrate as nitrogen source, remarl<ably the K. iactis 
KicarIA KLLA0F27995g/Kigbu1A mutant (IMK560) exhib- 
ited a 22.5-fold lower specific growth rate on guanidinobu- 
tyrate as sole nitrogen source (c. 0.02 h"^) (Fig. 6). 
Additionally, the strain IMK560 (KicarIA KLLA0F27995g/ 
KigbulA) was still able to grow on medium with arginine 
as sole nitrogen source at a specific growth rate of 
0.026 h"^ (versus 0.065 h"^ for IMS0367 in shal<e flask 
culture). In line with the labelling experiment, these results 
reinforced the predominant role of a pathway that involves 
an early transamination and the mobilization of the pep- 
tidic amino group which can still supply nitrogen for 
growth. The reduced specific growth rate of the strain 
IMS0367 on arginine (0.065 h"^) relative to growth on 
guanidinobutyrate (0.45 h"^) also suggests that the flux 
through the alternative arginine degradation pathway is 



limited by one or more reactions upstream the guanidi- 
nobutyrase (arginine transamination, ketoarginine decar- 
boxylase and/or oxidation of guanidinobutanal) (Fig. 5). 

Phylogenetic analysis of ureohydroiases 

While guanidinobutyrase activity had previously been 
measured in Aspergillus niger (Kumar and Punekar, 
2014), no gene encoding this activity has been identified 
in eukaryotes, a phylogenetic analysis was performed by 
searching for KICarl (arginase) and KIGbul (guanidi- 
nobutyrase) orthologues among the predicted protein 
sequences encoded in the genomes of yeasts belonging 
to the Saccharomycotina group, a subphylum of the 
Ascomycota (Kurtzman, 2003; Dujon, 2010). Using a 
phylogenetic tree topology proposed by Dujon (2010) for 
the Saccharomycotina group, orthologues of KICarl 
were found in all interrogated yeast proteomes. In con- 
trast, KIGbul (guanidinobutyrase) orthologues were not 
found in Saccharomytacea genera such as Saccharomy- 
ces sensu stricto, Kazaciistania, Naumovozyma, Naka- 
seomyces and Vanderwaltozyma, proposed to have 
arisen after the WGD event that played a major role 
in the evolution of this subphylum. Remarkably, one 
out of the 34 annotated S. cerevisiae genome 
sequences, available from the SGD database (http:// 
www.yeastgenome.org/) belonging to strain ECU 18, did 
harbour a gene whose predicted protein sequence 
shared 62% similarity with KIGbul (Fig. 7). A similar 
sequence occurs in the S. cerevisiae wine strain VIN7. 
Further inspection of the genome structure of the wine 
strain EC1118 revealed that the KIGbul orthologue was 
found on a genomic region acquired by horizontal trans- 
fer, most probably from Torulaspora delbrueckil. This 
gene origin fits the phylogenetic distribution of KIGbul 
orthologues in Saccharomycatae genera (Zygosaccharo- 
myces, Lachancea, Torulaspora, Kluyveromyces and 
Eremotheclum), and in genera belonging to the CTG 
group and Dipodascacae (Kurtzman, 2003; Dujon, 2010) 
(Fig. 7). 

To broaden the search for possible fungal orthologous 
of KIGbul, we explored sequence data from the fungal 
phyla Basidiomycota, Ascomycota (excluding the Saccha- 
romycotina already presented above) and Zygomycota 
which, together, comprise more than 85% of the charac- 
terized fungi. Blastp analysis using the KIGbul amino acid 
sequence as query confirmed that the presence of 
orthologous genes was widespread among fungi. This 
screen was based on stringent sequence similarity criteria 
(Blastp score > 250 and E-value < 1 .OE-70); as a refer- 
ence, the KICarl amino acid sequence queried with 
KIGbul sequence exhibited Blastp scores of 58 and an 
E-value of 8.0E-10, much lower than for the KIGbul ortho- 
logues identified in species reported in Table S4. These 
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results suggest that guanidinobutyrase and, possibly, the 
newly identified alternative arginine degradation pathway 
are not limited to K. lactis but may be widespread among 
fungi. 

Discussion 

Characterization of an arginase independent arginine 
degradation pathway 

The identification of an alternative, arginase-independent 
pathway for arginine metabolism in K. lactis could 
only be achieved by a systems biology approach that 



involved inspecting and integrating data from various 
cellular levels (Mustacchi etal., 2006; Kim etal., 2012). 
Together these approaches provided a strong basis for 
proposing a new fungal pathway for arginine degrada- 
tion (Fig. 8). This pathway starts with a transamination of 
arginine into ketoarginine. Biochemical characterization 
of this enzymatic step was limited by the fact that 
ketoarginine is not commercially available and by the 
promiscuous substrate specificity of the numerous 
transaminases in K. lactis (Dujon etal., 2004). However, 
detection of a specific labelling pattern corresponding to 
^^Ce^^Na (M+^) arginine in the labelling experiments is 
fully consistent with removal and reincorporation of a 
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Fig. 8. Updated arginine catabolism in K. lactis. Thick lines indicate the new guanidinobutyrase pathway linking arginine to succinate. The 
grey lines represent the reactions absent in K. lactis but present in other eukaryotes (as shown in Fig. 1). Reactions: EC 3.5.3.1: arginase, EC 
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labelled/unlabelled primary amino group via a reversible 
transamination reaction. This arginine isotopologue was 
observed in strain K. lactis ll\/IS0367 {KIcarIA Arg+) as 
well as in the reference strain GG1632, which suggested 
that the reaction involved one or more constitutively 
expressed transaminases. In the second and third steps, 
ketoarginine is decarboxylated to guanidinobutanal and 
oxidized to guanidinobutyrate respectively. This coin- 
cided with the moderate transcriptional upregulation, 
(+2. 7-fold) of KLLA0E02707g, a gene that is highly 
similar to S. cerevisiae ARO10, which encodes a broad 
substrate TPP-dependent 2-oxo acid decarboxylase 
(Vuralhan etal., 2005; Romagnoli etal., 2012; Bolat 
etal., 2013), a key component of the Erhlich pathway 
(Hazelwood etal., 2008). There are no literature data on 
activity of fungal TPP-dependent 2-oxo acid decarboxy- 
lases with ketoarginine as a substrate in yeasts. 
However, ketoarginine decarboxylase activity has been 



demonstrated in Pseudomonas species (Stalon and 
Mercenier, 1984; Yang and Lu, 2007b). The next step 
corresponds to a ureohydrolase reaction that converts 
guanidinobutyrate into y-aminobutyrate and participates 
in the mobilization of nitrogen from the arginine skeleton 
by releasing urea. The role of the KIGBU1 encoded 
protein has been demonstrated both in vitro and in vivo. 
Finally, y-aminobutyrate can be converted into succinate 
by the sequential activity of a GABA transaminase and 
succinate semialdehyde dehydrogenase. A transcript 
corresponding to KLLAF20548g gene, which shows 
strong sequence similarity with the S. cerevisiae GABA 
transaminase, was significantly upregulated in K. lactis 
IMS0367 (+4.7-fold). The topology of this proposed 
arginine catabolic pathway in K. lactis is similar to the 
arginine:pyruvate transaminase pathway, one of the 
several arginine degradation routes described in Pseu- 
domonas aeruginosa (Yang and Lu, 2007a,b). 
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Guanidinobutyrase is widely distributed in Fungi 

Most of the components of the proposed new pathway for 
arginine catabolism in K. iactis are also involved in other 
parts of metabolism: the transaminases are primarily 
involved in amino acid biosynthesis, the 2-oxo acid decar- 
boxylases are involved in degradation of certain amino 
acids (Hazelwood etal., 2008), while y-aminobutyrate 
transaminase and succinate semialdehyde dehydroge- 
nase are involved in metabolism of GABA, which can be 
formed independently of guanidinobutyrate through decar- 
boxylation of glutamate (Bach etai., 2009). The only 
enzyme activity that appears to be specific to this alterna- 
tive arginine degradation pathway is the guanidinobu- 
tyrase. We identified KIGbul orthologues in a wide 
diversity of fungal genomes. Together, these sequences 
form a distinct group of putative guanidino-acid ureohydro- 
lase genes that is well separated from genes encoding 
guanidino-amino acid hydrolase (arginase) and guanidino- 
amide (agmatinase) ureohydrolase groups. Several 
members (Table S4) of the new guanidino-acid ureohydro- 
lase group were previously annotated as agmatinase [e.g. 
Sclieffersomyces {Pichia) stipitis {SPEB2) agmatine ureo- 
hydrolase XM_001 385297]. The present study indicates 
that functional annotation of these genes may have to be 
revised (Jeffries etal., 2007; Smith etal., 2008). Heterolo- 
gous expression in a carl A strain of S. cerevisiae, followed 
by growth tests on media with either arginine or guanidi- 
nobutyrate as sole nitrogen source, offers a relatively 
simple approach to investigate whether these genes 
indeed encode guanidinobutyrases. 

To further explore the diversity of fungal arginine catabo- 
lism, we investigated whether arginine decarboxylase 
genes could be identified in fungi. This would enable an 
alternative to the pathway described above, since arginine 
decarboxylase yields agmatine, the substrate for agmati- 
nase, whose decarboxylation yields putrescine. This 
pathway is found in bacteria (Szumanski and Boyle, 1990; 
Sekowska etal., 2000), plants (Klein etal., 1999) and 
animals (Iyer etal., 2002). A Blastp search using the 
Escherichia coli arginine decarboxylase SpeA amino 
acid sequence as a query failed to identify predicted 
protein sequences in the fungal kingdom (thresholds 
E-value < 1.0E-20, Blast score > 100). We could hypoth- 
esize that the occurrence of a guanidino-acid ureohydro- 
lase is correlated with the absence of arginine 
decarboxylase and consequently of a guanidino-amide 
ureohydrolase. This would indicate that fungi synthesize 
putrescine directly by decarboxylation of ornithine exclu- 
sively. Therefore the presence of guanidino acid ureohy- 
drolase pathway would then represent an alternative to 
arginase in absence of an agmatinase pathway. 

This new pathway described in this study, links arginine 
to succinate and demonstrates different entry points in 



central carbon and nitrogen metabolisms. It definitely 
increased the nitrogen source spectrum of K. lactis ena- 
bling this yeast to grow on medium containing guanidi- 
nobutyrate as sole nitrogen source; however, the 
physiological significance of this pathway remains unclear. 
Dedicated studies to investigate transcriptional regulation, 
localization of the enzymes comprising the pathway and 
connection to the arginine metabolism (arginase pathway 
and the polyamine synthesis) are necessary to understand 
the metabolic, regulatory and evolutionary meanings of the 
guanidinobutyrase pathway deeper. 

Guanidinobutyrase has been lost after the WGD event 

Although widespread in fungi, guanidinobutyrase genes 
were absent in Saccharomycetaceae clades (Saccharo- 
myces sensu stricto, Vanderwaltozyma, Nakasoemyces, 
Naumovozyma) (Fig. 7) that arose after the WGD, an 
iconic event in the recent evolutionary history of the Sac- 
charomycetaceae that occurred nearly 100 million year 
ago (Wolfe and Shields, 1997; Seoighe and Wolfe, 1999). 
The WGD is proposed to have involved the transient 
formation of a tetraploid ancestor, which gradually 
returned to a normal ploidy (2n) through a combination of 
a large number of deletion and epigenetic silencing 
(Sankoff etal., 2012; Zheng and Sankoff, 2012). The 
current model proposes that S. cerevisiae originates from 
an ancestor with a genome of ~ 10 000 genes, a number 
that has been reduced to less than 6000 genes in S. 
cerevisiae, of which 1104 form 552 so-called ohnolog 
pairs (Scannell etal., 2006; Byrne and Wolfe, 2007). The 
guanidinobutyrase example is consistent with the notion 
that, for some genes, both ohnologs may have been lost 
during post-WGD genome size reduction. The observa- 
tion that this particular reaction has been lost during 
genome size reduction may, after further physiological 
analysis of the enzymology and regulation of different 
fungal pathways for arginine catabolism, provide clues on 
the environmental conditions under which this event took 
place. 

Interestingly, a few S. cerews/ae strains appear to have 
re-acquired KIGbul orthologue (Novo etal., 2009). 
Recent genome sequence analysis of the wine strains 
EC1118 (Novo etal., 2009) and VIN7 (Borneman etal., 
2012) demonstrated that these strains have, by horizontal 
gene transfer, acquired DNA segments from microorgan- 
isms that share their natural habitats. In both strains, the 
KIGBU1 orthologous genes were found in an introgres- 
sion region (region A, 38 kb) on chromosome VI (Novo 
etal., 2009) that most probably originated from T. del- 
brueckii (Gordon etai, 2011). Arginine is together with 
proline, the most abundant nitrogenous compound in 
grape must (Beltran etal., 2005). Under anaerobic condi- 
tions, arginine catabolism via the arginase pathway 
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Table 6. Kluyveromyces lactis and Saccharomyces cerevisiae strains used in this study. 



Strain 


Organism 


Genotype 


Reference 




K. 


iactis 


Prototrophic reterence strain 


uujon Gi at. yd.uv^) 


OO 1 DO^ 


K. 


iactis 


l\IKUoU/\..IOXr 


r\,OOIStra Gi 31. [eiyJU^) 


MVi r-J\od. 


K. 


iactis 


f\lKuoUl\..IOXr f\LLMUrUU/U^A \i\lCar t i\j..lOXr-KcinivtA-tOXr 


This study 


MVIr\4oo 


K. 


iactis 


l\IKUoU/\..IOXr r\LLAUDUO£i4/A [t\ICart::A}..IOXr-KaniVIA-IOXr 


This study 


MVI rx'fo'l- 


K. 


iactis 


i\lt\UOUiX..lOXr fVLLMUU lU/oDZi \f\ipr00iA} ..lOXr-KanivlA-IOXr 


This study 


MvloUob/ 


K. 


iactis 


l\IKUoU/\..IOXr r\LLAUrUU/U4A [l\ICari/\)..IOXr-KaniVIA-IOXr f\\Q 


This study 


MvirxDDU 


K. 


iactis 


r\lt\UOUl\..IUXr rxLLMUruU/ UtA llAj..llJXr-nctlllVIA-IUXr 


1 [lib biuuy 








i\lyDU tOXr-amOS Y IVI-IOXr 




pcM pk'1 1 7 n 


S. 


cerevisiae 


D frttrttKrtrih 1/^ rof Q i*(!iri/^Q fl^dTo 

r ruiuirupiiic [yiyryriuy bLidiii iviMia 




r^PM DLCi-lo c: n 
OtlN.rtxl 1 O-O U 


S. 


cerevisiae 


MM la uraJ-O^ 


tn+ian anH l^nHar IOC\f\'7\ 

cniian anu r\oiiBr 


MVI rv'fOO 


S. 


cerevisiae 


ivIMIa UraJ-O^ Cat 1 /l..lOXr-KaniVlA-IOXr 


This study 


MVI rv'foO 


S. 


cerevisiae 


ivIMia Urao-Oi^ Caf£LlA..lOXr-KanMA-lOXr 


This study 


l^/lU'/1/1 1; 

MVI rv'f^O 


S. 


cerevisiae 


iviMia ursLo-Oi^ prooiA..ioxr-KaniviA-toxr 


This study 


IIVIZSIO 


S. 


cerevisiae 


MATa ura3-52 carmr.loxP-lnanMX-ioxP pUDE168 


This study 








(TDH3pr.KLCAR1-CYC1,er URA3 2n) 




IIVIZSH 


S. 


cerevisiae 


MATa ura3-52 car2A.;toxP-/canMX-/oxP pUDE169 


This study 








(TDH3p,.KLCAR2-CYC1,er URA3 2^) 




IMZ342 


S. 


cerevisiae 


MATa ura3-52 pro3A::ioxP-l<ani\/IX-ioxP pUDE170 


This study 








(TDH3pr.KLPR03-CYC1,er URA3 2^) 




IMZ312 


S. 


cerevisiae 


MATa ura3-52 pAG426GPD-ccdB {TDH3prCYC1 URA3 2|i) 


This study 


IIV1E215 


S. 


cerevisiae 


MATa ura3-52 pUDE264 {TDH3prKIGBU1-CYC1,e, URA3 2ix) 


This study 


ll\/IE216 


S. 


cerevisiae 


MATa ura3-52 PUDE265 {TDH3prKLLA0F04235g-CYC1,„ URA3 2\i) 


This study 



enables utilization of tliree of its four nitrogen atoms, two 
as a result of the arginase reaction and one from the 
ornithine transaminase step (Figs 1 and 8). The fourth 
nitrogen atom ends up in proline which, under anaerobic 
conditions cannot be assimilated by S. cerevisiae, as 
proline catabolism requires molecular oxygen. In contrast, 
arginine catabolism through the guanidinobutyrase 
pathway should enable the mobilization of all four nitrogen 
atoms, one issued from the arginine transamination, two 
from the guanidinobutyrase reaction, and the last one 
from GABA transamination. Therefore, the 'guanidinobu- 
tyrase pathway' could provide a selective advantage 
under the anaerobic and nitrogen limited conditions of 
wine fermentation (Bisson, 1999). The present study 
therefore provides a clear lead for exploring the signifi- 
cance of a horizontal gene transfer event in wine yeasts 
that can be explored by functional analysis of the wine- 
yeast KIGBU1 orthologue. 

The present study clearly demonstrates that in K. lactis 
and more generally in Fungi arginine degradation was not 
limited to the ubiquitous arginase pathway. By using a 
top-down approach we could identify a new fungal meta- 
bolic pathway involved in arginine degradation implicating 
a guanidinobutyrase, an enzyme activity that could be 
widespread in Fungi. 

Experimental procedures 

Yeast strains and maintenance 

The S. cerevisiae and K. lactis strains used in this study are 
listed in Table 6. The S. cerevisiae strains were constructed in 



the CEN.PK background (Entian and Kotter, 2007; Nijkamp 
etal., 2012) and the K. /ac//s strains were all derived from the 
strain CBS2359 (Wesolowski-Louvel etal., 1996; Dujon 
etal., 2004; Kooistra etal., 2004). Yeast strains were main- 
tained on YPD medium (demineralized water; 10 g yeast 
extract; 20 g peptone; 20 g I ' glucose). Culture stocks 
were prepared from shake flask cultures incubated at 30°C 
and stirred at 200 r.p.m., by addition of 20% (v/v) glycerol and 
were stored at -SOX. 

Deletion of CAR 1, CAR2, PR03 and KLLA0F27995g 
(KIGBU1 ) in S. cerevisiae and K. lactis 

Gene deletions in S. cerevisiae were performed by Integra- 
tion of loxP-kan-loxP (kanMX) cassettes (Guldener etal., 
1996) via the short-flanking-homology PCR method (Wach 
etal., 1994). Sequences of oligonucleotide primers are 
shown in Table S5. Deletion cassettes for CAR1, CAR2 and 
PROS were amplified using Phusion Hot-Start polymerase 
(Finnzymes, Landsmeer, the Netherlands) and the template 
plasmid pUG6 (Guldener etal., 1996) using primers 
ScCARI Fw/ScCAR1 Rv, ScCAR2Fw/ScCAR2Rv and 
ScPROSFw/ScPROSRv respectively. The transformation of 
S. cerews/ae CEN.PK11 3-5 D with the CAR1, CAR2orPR03 
deletion cassettes were performed using the LiAc method as 
previously described in Gietz and Woods (2002) resulting in 
strains IMK435, IMK436 and IMK445 respectively. Cassettes 
for deletion of the corresponding K. lactis genes were simi- 
larly obtained from pUG6 using primers KICARIFw/ 
KICARIRv, KICAR2Fw/KICAR2Rv and KIPR03Fw/ 
KIPR03RV. The resulting deletion cassettes were 
transformed into K. iactis GG1632 resulting in strains 
IMK432, IMK433 and IMK434 respectively (Kooistra etal., 
2004). The deletion of the gene KLLA0F27995g (KIGBU1) 
was performed by the integration of a loxP-amdsYM-loxP 
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Table 7. Plasmids used in this study. 



Plasmid 


Characteristic 


Reference 




PGR template tor loxP-K3nMX4-loxP cassette 


(^iilHon/^i- at o/ /iOQC\ 

ciuiaener ei ai. ^lyyo^ 


nl in hnhNTI 


milipidLtj lUi luXr-lipfllv I / -/CJAr OdobcLlt; 


yjUG\Uti\\ii\ til at. ^^UU^^ 


nl in AmHcV 


rOr\ milipidLtj lUl /UAr-cI(//UO / /W-ZCJAf^ UdbbtJltc 




pon-zfcsru-uiuiii 






nUD180 


Gateway entry clone, \^CAH1 


This study 


pUD181 


Gateway entry clone, K\CAR2 


This study 


pUD186 


Gateway entry clone, K\PR03 


This study 


pAG426GPDccdB 


2\x ori \JRA3 TDH3p,-ccdB-CYC1,er 


Alberti ef a/. (2007) 


PUDE168 


2|i ori URA3 TDH3p,-KICAR1-CYC1,„ 


This study 


PUDE169 


21^ ori URA3 TDH3p,-KICAR2-CYC1,sr 


This study 


PUDE170 


2n ori URA3 TDH3p,-KIPR03-CYC1 ,„ 


This study 


pUDE264 


2\^ ori TDH3prKIGBU1-CYC1,er 


This study 


pUDE265 


2n ori TDHSpr KLLA0F04235-CYC1 ,„ 


This study 



(Solis-Escalante efa/., 2013) cassette by homologous 
recombination at the locus of the corresponding gene using a 
cassette with long flanking homology region. The initial cas- 
sette was amplified from plasmid pUGamdsYM plasmid 
(Solis-Escalante efa/., 2013) using primer KIGBUIdell and 
KIGBU1del2, while the two flanking regions were amplified 
directly from genomic DNA using primers KIGBUIdelS/ 
KIGBU1del3 for the 5' UTR and primers KIGBU1del4/ 
KIGBUIdel? for the 3' UTR. The three resulting PGR 
products were purified and used as template for a fusion PGR 
with the nested primers KIGBU1del6/KIGBU1del8. The frag- 
ment with the appropriate length was gel purified and used to 
transform the strain 1MS0367 resulting in strain IMZ449 
(Kooistra efa/., 2004). 

Gorrect integration of the kanMX or amdsYM cassettes 
and replacement of the gene of interest was, in all cases, 
verified by diagnostic PGR using two sets of primers: a 
forward primer specific for the 5' UTR (untranslated region) of 
the gene (Table S5 labelled with conFw) and two reverse 
primers, one specific for the original gene (labelled with 
conRv) and one for the deletion cassette (tagged with KanA 
or amdSY). 

Cloning and overexpression of KIGAR1, KICAR2, 
KIPR03, KLLA0F27995gMGB[j: and KLL0F04235g in 
S. cerevisiae 

Genomic DNA of the prototrophic reference strains S. cerevi- 
siae GEN.PK113-7D and K. lactis GBS2359 was pre- 
pared as described previously (Burke efa/., 2000). ORFs 
KLLA0F00704g (KICAR1), KLLA0B05247g (KICAR2), 
KLLA0D10736g (KIPR03) KLLA0F27995g (KIGBU1) 
and KLLA0F04235g were cloned from genomic DNA 
using Phusion Hot-Start polymerase (Finnzymes) and pri- 
mers KIGARIclFw/KIGARIclRv, KIGAR2clFw/KIGAR2clRv, 
KIPR03clFw/KIPR03clRv, KIGBUIclFw/KlclGBUI Rv and 
KLLA0F04235gclFw/KLLA0F04235gclRv respectively. 

The first three PGR products were cloned into the pGR- 
zero-blunt vector (Invitrogen, Breda, the Netherlands) result- 
ing in plasmid pUDISO, pUD181 and pUD186 respectively. 
These plasmids as well as the PGR product deriving from the 
KIGBU1 and KLLA0F04235g genes and the p426GPDccdB 
(Alberti efa/., 2007) plasmid were digested with Spel and 



Xhol, purified from gel and llgated Into the expression vectors 
to yield PUDE168, pUDE169, pUDE170, pUDE264 and 
pUDE265 respectively (Table 7). The pUDE168 plasmid was 
transformed in strain IMK435 resulting in IMZ310; the 
pUD169 in strain IMK436 resulting in strain IMZ311; plasmid 
pUD170 in strain IMK445 resulting in strain IMZ342 and 
finally plasmid pUDE264 and pUDE265 in strain GEN.PK113- 
5 D resulting in strain IME215 and IME216 respectively. 



Media and culture conditions 

Growth experiments were conducted in synthetic medium 
containing salts, trace elements and vitamins, prepared and 
sterilized as described previously (Verduyn efa/., 1992). 
Glucose was added to a final concentration of 20.0 g ^^ 
When ammonium sulphate was not the nitrogen source in the 
synthetic medium, it was replaced by L-arginine, L-ornithine, 
L-proline or L-guanidinobutyrate, these compounds were filter 
sterilized and added to sterile medium to concentrations of 
3.3 g I \ 5.0 g I \ 8.7 g f and 2.9 g respectively. Moreover, 
3.3 g f potassium sulphate was added to compensate for the 
removal of ammonium sulphate. If required, 0.15 g uracil 
and/or 200 mg of G418 were added to complete media. 
Selection Agar plates were made by adding 20.0 g I ' agar to 
these synthetic media. Agar plates for selection of the 
amdsYM marker were prepared as described previously 
(Solis-Escalante etal., 2013). 

Shake flask cultures were conducted In 500 ml or 100 ml 
shake flasks containing 100 ml or 20 ml of liquid medium 
respectively and incubated in an orbital shaker (New Brun- 
swick Scientific, Edison, NJ) at 200 r.p.m. at 30°G. Gontrolled 
aerobic batch cultures were grown in 2 I laboratory bioreac- 
tors (Applikon, Schiedam, the Netherlands) with a working 
volume of 1.2 I. Gulture pH was maintained at pH 6.0 by 
automatic addition of 2 M KOH (ADI 1030 Biocontroller, App- 
likon) and temperature was kept constant at 30°G (Lauda 
E300 thermostat, Lauda DR. R. Wosber Gmbh & Go. KG 
Lauda-Konigshofen, Germany). Aeration was achieved by 
sparging air through the cultures at a flow rate of 0.5 I min"' 
continuously stirring at 800 r.p.m. The exhaust gas of the 
bioreactor was cooled in a condenser using a cryostat set at 
4°G (Lauda El 00 cryostat). GO2 and O2 concentrations In the 
off gas were measured with an NGA 2000 analyser (Rose- 
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mount Analytical, Solon, OH). For spot plate assays, pre- 
cultures in synthetic medium with glucose diluted to 10*^ cells 
m\ \ which corresponded to an ODeeo of 0.1 (Hazel wood 
etal., 2006). Serial dilutions (10''-10^ cells mM) were made 
and 5 |xl of each dilution was spotted on agar plates and 
incubated at 30°C for 3 days. 



Preparation of cell extracts 

For preparation of cell extracts, culture samples were har- 
vested by centrifugation, washed twice with 10 mM potas- 
sium phosphate buffer (pH 7.5) containing 2 mM EDTA and 
stored at -20°C. Before cell disruption, samples were thawed 
at room temperature, washed, and resuspended in 100 mM 
potassium phosphate buffer (pH 7.5) containing 2 mM MgCl2 
and 2 mM dithiothreitol. Extracts were prepared by sonication 
with 0.7 mm glass beads at 0°C for 2 min at 0.5 min intervals 
with an MSB sonicator (Wolf Laboratories Limited, Pockling- 
ton, UK) (150 W output; 8 |xm peak-to-peak amplitude) (Luttik 
etal., 2008). Unbroken cells and debris were removed by 
centrifugation at 4°C (20 min; 36 000 g). The resulting cell 
extract was used for enzyme assays. 



Enzyme activity assays 

For the arginase enzymatic assay 50 ^1 cell extract were 
activated in 950 nl manganese maleate buffer (50 mM man- 
ganese sulphate, 50 mM maleic acid, pH 7) for 1 h at 37°C 
(Messenguy etal., 1971). The reaction mixture for arginase 
assays, prepared in dark Eppendorf tubes, contained 60 |xl of 
activated cell extract, 400 ^1 of 713 mM arginine solution (pH 
9.5) and demineralized water up to 1 ml. The reaction mixture 
was incubated for 30 min at 37°C. To stop the reaction, 0.7 ml 
sulphuric-phosphoric acid mixture (20% v/v concentrated sul- 
phuric acid and 60% v/v syrupy phosphoric acid in demineral- 
ized water) was added to the reaction mixture. The amount 
of urea produced was quantified using the Archibald 
method (Archibald, 1945) with a calibration line ranging 
from 0 until 0.6 mM of urea. 0.06 ml of a 4% v/v 
a-isonitrosopropiophenone in ethanol solution was added and 
samples were thoroughly mixed before boiling for 1 h in a 
100°C water bath to develop the colour. The samples were 
cooled at room temperature for 15 min and the absorbance at 
540 nm was measured in a Libra S11 spectrophotometer 
(Biochrom, Cambridge, UK). 

The reaction mixture for ornithine amino transferase (OAT) 
was prepared in dark eppendorf tubes, containing, in a 1 ml 
final volume: 50 mM ornithine, 20 mM a-ketoglutarate, 1 mM 
pyridoxal 5-phosphate and 100 mM potassium phosphate 
buffer pH 8 (Vogel and Kopac, 1960) and 2.5 to 50 |xl cell 
extract. The mixture was incubated for 30 min at 37°C before 
the reaction was stopped by addition of 0.5 ml 10% trichlo- 
roacetic acid. The amount of pyrroline-5-carboxylate pro- 
duced was determined by the O-amino-benzaldehyde assay 
using a calibration line ranging from 0 to 0.5 mM. Colour was 
developed by addition of 0.5 ml 0.5% o-aminobenzaldehyde 
in 95% ethanol and incubation at room temperature for 1 h. 
Samples were then centrifuged and the clear supernatant 
was used to measure absorbance at 440 nm with a Libra S11 
spectrophotometer (Biochrom). 



The reaction mixture for pyrroline-5-carboxylate (P5C) 
reductase contained 50 mM A/-[Tris(hydroxymethyl)methyl]- 
3-aminopropanesulphonic acid (TAPS), 0.14 mM NADH and 
cell extract in a final volume of 300 ^1 (Brandriss and 
Magasanik, 1979). Enzyme activity was monitored as NADH 
consumption at a wavelength of 340 nm using a GENios 
Microplate Reader (Tecan, Giessen, the Netherlands), after 
addition of pyrroline-5-carboxylate at pH 7 to a final concen- 
tration of 1 mM. A^-pyrroline-5-carboxylate (P5C) was syn- 
thesized through oxidation of 5-hydroxylysine as previously 
described in Williams and Frank (1975). To start, 71.4 mM 
hydroxylysine was prepared in a brown bottle. 50 mM sodium 
metaperiodate was prepared in dim light, neutralized with 1 M 
NaOH and 4.4 ml of the neutralized pre-cooled periodate 
solution were quickly added to 2.8 ml of ice-cold hydroxyly- 
sine solution and incubated for 8 min at 4°C. Remaining 
periodate was destroyed by addition of 70 |al of 1 M glycerol 
and incubation for 2 min. The mixture was then acidified by 
the addition of 60 |il of a 6 M HCI solution. The complete 
reaction mixture was applied to a 5 ml Dowex 50 column, 
which was equilibrated with water after thorough washing 
with a 1 M HCI solution. Elution of the P5C was performed by 
applying 100 ml of 1 M HCI and collecting 2 ml fractions. For 
qualitative detection of P5C 20 |xl of each elution fraction 
were mixed with 20 ^1 of 3 M sodium acetate and 500 |il 
0.15% ninhydrin in glacial acetic acid. The mixture was incu- 
bated on a hot plate at 100°C. If more than 0.5 ^g of P5C 
were present in the test tube the mixture turned cherry-red 
shortly before reaching the boiling point of the solution. 
Elution fractions showing presence of sufficient amounts of 
P5C were then analysed quantitatively. Aliquots of the elution 
fractions containing P5C were mixed with HCI to a final con- 
centration of 1 M HCI. Subsequently 200 ^1 of this mixture 
were added to 100 |xl of 0.5% o-aminobenzaldehyde solution 
in 100% ethanol. After incubation for 40 min at room tempera- 
ture the absorbance was measured at 444 nm. Concentra- 
tions were calculated using the molar absorption coefficient 
of 2940 M-^ cm ^ 

The reaction mixture for guanidinobutyrase (GBU) 
enzyme assays was prepared in dark eppendorf tubes, con- 
taining in a 1 ml final volume: 50 mM glycine buffer (pH 9), 
5 mM MnS04 and 50 ^1 to 100 nl cell extract. The reaction 
was started by addition of 50 mM guanidinobutyric acid. 
After 30 min of incubation at 37°C, the reaction was 
stopped by addition of 700 |xl of sulphuric-phosphoric acid 
mixture (20% v/v concentrated sulphuric acid and 60% v/v 
syrupy phosphoric acid in demineralized water). The 
amount of urea produced was quantified using the Archibald 
method (Archibald, 1945) with a calibration line ranging 
from 0 until 0.6 mM of urea. 0.06 ml of a 4% v/v 
a-isonitrosopropiophenone in ethanol solution was added 
and samples were thoroughly mixed before boiling for 1 h in 
a 100°C water bath to develop the colour. The samples 
were cooled at room temperature for 15 min and the 
absorbance at 540 nm was measured with a Libra S11 
spectrophotometer (Biochrom). 

Isotopomer labelling 

The biomass was pre-cultivated in 500 ml shake flasks con- 
taining 120 ml of glucose synthetic medium with 3.3 g I"' of 
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non-labelled arginine as sole N-source. During the exponen- 
tial growth phase the cells were harvested, washed with 
10 mM potassium phosphate buffer (pH 7.5) and subse- 
quently suspended in 2 separate shake flasks containing 
50 ml minimal medium without any nitrogen source. Immedi- 
ately after the addition of the labelled and unlabelled arginine 
to the 2 shake flasks, an initial sample of 800 |il was taken for 
intracellular metabolite measurement as described below. 
Samples were taken over the course of the experiment for 
half cell-duplication time. Intracellular metabolite samples 
obtained from cultivations without labelled material were 
used to determine intracellular metabolite concentrations. 
Isotopomer measurement data derived from the shake flasks 
containing the ^^C6'^N4 labelled arginine was used for meta- 
bolic flux analysis. 

Sampling and analytical techniques 

Specific growth rates were estimated from the increase in 
optical density at 660 nm (ODeeo), as measured with a 
Novaspec II spectrophotometer (GE Healthcare Europe 
GmbH, Diegem, Belgium) after dilution of culture samples to 
ODeeo values below 0.3. Biomass dry weight was determined 
by filtrating 10 ml culture samples over dry, pre-weighed nitro- 
cellulose filters (Gelman laboratory, Ann Arbor, Ml). After 
washing with demineralized water, filters were dried using a 
360 W microwave oven and immediately weighed. Using 
multiple OD and corresponding dry weight measurements, a 
stoichiometric relation between these two parameters was 
determined for each of the batch cultures. 

Extracellular metabolite samples were obtained after cen- 
trifugation of the culture broth. The supernatant was subse- 
quently analysed using the analytical technology described 
below. 

Glucose and ethanol concentration were determined using 
HPLG (Waters Alliance 2695 Separation Module combined 
with a Waters 2487 Dual Absorbance Detector and a Waters 
2410 Refractive Index Detector, Waters, Milford, MA). Elution 
was achieved using a Bio Rad HPX87H column (Bio-Rad) at 
60°C with 0.5 mM H2SO4 at a flow rate of 0.6 ml-min ^ Free 
amino acid concentrations were determined using the AccQ- 
Tag Ultra Derivatization Kit (Waters) and measurement by 
HPLG (Waters Alliance 2695 Separation Module combined 
with a Waters 474 fluorescent detector. Waters) according to 
the procedure described in vanWandelen and Gohen (1997). 
Extracellular urea and ammonia was measured using an 
enzyme essay kit (K-URAMR, Megazyme, Ireland). Using 
this kit the amount of both compounds was quantified by 
measuring the decrease in absorbance at 340 nm in a 
GENios Microplate Reader (Tecan). 

Intracellular metabolite sampling requires rapid quenching 
of the cells using -40°C 80%/20% methanol/ethanol mixture. 
After extraction the 800 ^il samples were added directly to 5 ml 
tubes of pre-cooled methanol. Subsequently cold methanol 
was used to wash the samples and the metabolites were 
extracted with boiling ethanol as previously described by 
Ganelas etal. (2008). Accurate intracellular metabolite quan- 
tification required the addition of an internal standard in the 
form of uniformly labelled ^^G-cell-extract (Wu etal., 2005). 
Using GG-MS, the concentration and mass shift of ^^C^U 
samples was measured for, pyruvate (Pyr), alanine (Ala), 



glycine (Gly), valine (Val), leucine (Leu), iso-leucine (lie), 
proline (Pro), serine (Ser), threonine (Thr), methionine (Met), 
aspartate (Asp), phenylalanine (Phe), cysteine (Gys), gluta- 
mate (Glu), lysine (Lys), asparginine (Asp), glutamine (Gin), 
tyrosine (Tyr), histidine (His), ornithine (Orn) and tryptophan 
(Trp) using the protocol described by de Jonge etal. (2011). 
GC-MS was also used to quantify iso-citrate (i-Cit), citrate 
(Git), malate (Mai), fumarate (Fum), a-ketoglutarate (aKg), 
putrescine (Put) and y-aminobutyrate (GABA) concentrations 
and '^C'^N enrichment (Cipollina et al., 2009). Succinate (Sue) 
was measured using LC-MS according to methods previously 
developed by van Dam etal. (2002). Arginine and ornithine 
was measured by ion pair reversed phase liquid chromatog- 
raphy combined with tandem mass spectrometry (Seifar et al., 
2009) after derivatization using o-phetalaldehyde/isobutyryl- 
L-cysteine (OPA-IBLC) in accordance to Hess (2012). The 
data derived from the isotopomer measurements were cor- 
rected for naturally occurring isotopes as well as mass-shifts 
caused by derivatization and non-carbon isotopes. 



"C'^A/ metabolic flux analysis 

All calculations were done using the gPROMS software 
suite (Process Systems Enterprise Limited, London, UK) 
in which all fluxes and intracellular concentration were 
assumed to be in steady state. The set of differential equations 
consisted of two distinct parts; initially the fate of the ^^N in 
nitrogen containing compounds was determined assuming a 
sole influx of labelled material through arginine. For each 
measured amino acid the corresponding incorporation into 
new biomass was calculated assuming the biomass compo- 
sition as reported by Duarte et al. (2004). The '^C'^^N metabolic 
model established in this study contains 30 metabolites and 46 
metabolic reactions. A more elaborate description of model 
parameters, data fitting and assumptions can be found in the 
supplemental data (Tables S2 and S3; Figs SI and S2). 



RNA-seq transcriptome analysis 

Biomass samples for RNA extractions were sampled from 
mid-exponential phase of K. lactis strains GG and IMS0367 
grown in batch cultivations in bioreactors with chemically 
defined medium with glucose and arginine as sole N-source. 
The sampling, the extraction and the cDNA synthesis was 
performed as previously described in Hazelwood et al. (201 0). 

Sequencing was performed using lllumina Genome Ana- 
lyser II and carried out by Baseclear (Leiden, the Netherlands). 
Data sets of 50 bp single reads of at least 1-Gb were gener- 
ated. RNA-seq data generated in this study were submitted to 
Genome Expression Omnibus database (http://www.ncbi 
.nlm.nih.gov/geo/) under Accession No. GSE56060. 

The genome sequence of K. lactis strain CBS2359 and its 
annotations were retrieved from the genolevures database 
(Martin etal., 2011) and used for all analysis. The data were 
aligned to the reference using Burrows-Wheeler Alignment 
(BWA) tool (Li and Durbin, 2009) and the resulting Sequence 
Alignment/Map (SAM) file was converted to a sorted and 
indexed Binary Alignment/Map (BAM) format using samtools 
(Li and Durbin, 2009; Li etal., 2009). Gene expression levels 
were estimated using FPKM values by the Cufflinks software 
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(Trapnell et al., 2010). To identify differential gene expression 
between strains GG1632 and IMS0367 grown in bioreactor 
batcli cultures with arginine as sole nitrogen source, RNAseq 
data comparison was performed and statistically assessed 
using Cuffdiff (Trapnell etal., 2010). 

Phylogenetic distribution anaiysis 

The phylogenetic tree of the Saccharomycotina was based 
on the previously postulated structure in Kurtzman (2003) 
and Dujon (2010). The sequence of KIGBU1 orthologous 
genes were identified by Blastp and tBIastn in the genome 
sequences of the microorganisms mentioned in Fig. 8. 
The Blast significance thresholds were set at an 
E-value < 1 .OE-70 and a Blast score > 250. 
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